Spin readout and initialization in a semiconductor quantum dot 
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Electron spin qubits in semiconductors are attractive from the viewpoint of long coherence times. 
However, single spin measurement is challenging. Several promising schemes incorporate ancillary 
tunnel couplings that may provide unwanted channels for decoherence. Here, we propose a novel 
spin-charge transduction scheme, converting spin information to orbital information within a single 
quantum dot by microwave excitation. The same quantum dot can be used for rapid initialization, 
gating, and readout. We present detailed modeling of such a device in silicon to confirm its feasibility. 
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0^ 



X3 

o 

> 

O 
m 
o 



I 

o 
o 



X 



Quantum dot qubits lie at the point of intersection be- 
tween quantum computing and few-electron spintronics. 
The long coherence times and novel quantum mechanical 
properties of electron spins, the attraction of scalability, 
and the shared technology base with the digital electron- 
ics industry fuel a large research effort to develop these 
systems. Several proposals have been put forward for 
spin qubits in semiconductors 0,|3j0>Sl3- In particular, 
exchange coupled gates in semiconductor dots appear to 
be workable |a| , and several components of qubit technol- 
ogy have been demonstrated recently However, the 
combined challenge of preparing, storing and measuring 
spins is formidable. 

Measurements of spin qubits pose a special challenge. 
On the one hand, qubits should be well isolated from 
their environment to avoid decoherence. On the other 
hand, it is necessary to individually couple the qubits 
to an external measurement device. Qubit initializa- 
tion involves an additional dissipative coupling to the 
environment. For quantum computing, we must initiate 
such couplings selectively, and with sufficient strength to 
perform the operations quickly. Indeed, scalable quan- 
tum computing relies on fault-tolerant quantum error 
correction algorithms, involving frequent, parallel mea- 
surements, and a steady supply of initialized qubits 
Fortunately, rapid and sensitive quantum measurement 
techniques involving radio frequency single electron tran- 
sistors (rf-SETs) have been developed [g- rf-SETs have 
been used to detect the tunneling of individual electrons 
in semiconductor devices . 

We have previously designed a quantum dot architec- 
ture specifically for the purpose of manipulating electron 
spins for fast and accurate two-qubit operations that 
serve as universal gates for quantum computations j^. 
Recent experimental results have shown that decoherence 
does n ot p ose a fundamental problem for such gate oper- 
ations . Using special qubit geometries , it should 
be possible to perform reliable gate operations in silicon 
quantum dots at rates between about 1 MHz and 1 GHz 
|l2l |. The question is now whether similar speeds and 
reliable operation can be achieved for measurement and 
initialization operations. In this paper, we show how to 



extend our architecture to enable the readout and initial- 
ization of the state of a single spin. The general scheme 
does not assume a specific material system. However, in 
this Letter we focus upon silicon-based devices because 
of their desirable coherence properties. 

Our method relies on the idea of converting spin in- 
formation to charge information discussed by Loss and 
DiVincenzo 1], and using single electron transistors to 
read out the resulting spin state as proposed by Kane et 
al. 2] , who posit spin-dependent charge motion onto im- 
purities in Si. Martin et al. have proposed a scheme 
for single spin readout that also converts spin informa- 
tion into charge information in an electron trap near a 
conducting channel. The resistance of the channel de- 
pends on the occupation of the trap, which in turn can 
be made to depend on the spin. Our proposed device 
has the additional feature in that it enables both readout 
and rapid initialization of the spin state. Rapid initial- 
ization (as opposed to initialization by thermalization) 
is just as essential as readout, and just as difficult to 
carry out. Our design is simple in that it obviates the 
need for spin-polarized leads or ancillary qubits. We 
go beyond conceptual design-the full three-dimensional 
system is simulated self-consistently and the operating 
parameters are thereby optimized. This enables us to 
demonstrate quantitatively that the device is in a work- 
ing regime. 

The basic Si-SiGe heterostructure is described in Ref. 
5]. The main point for present purposes is that the 
active layer is pure strained Si, which minimizes deco- 
herence from spin-phonon coupling jl4l |. The qubits are 
gated quantum dots which hold one electron, and the 
gate geometry is the key to our scheme. Specifically, 
the electrons are confined in asymmetric lateral wells, 
such that orbital excitation results in lateral center-of- 
charge movement. Figure 1 shows the first two orbital 
states of an electron confined to an asymmetric box, with 
center-of-charge positions varying by a distance Ay. If 
an external radiation source of frequency Ei2/h drives 
the system between the ground and first excited states, 
n = 1 and n = 2, the center-of-charge of the electron will 
oscillate in time at the Rabi frequency Vf;. This motion 
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FIG. 1: Schematic example of a one-dimensional, asymmetric 
confining potential. Quantum dot energy states and transi- 
tion rates for readout and initialization are shown. 



can be detected by a sensitive electrometer, such as a sin- 
gle electron transistor or a quantum point contact. With 
the system placed in a magnetic field, charge motion can 
be generated in a spin- dependent fashion. We define the 
logical and 1 of the qubit as the spin up and spin down 
states of the electron in its orbital ground state, and use 
the orbital excited state only during initialization and 
readout. The system is now driven at the readout fre- 
quency vi2{B) — Exijh — giJ,BB/h. This causes charge 
motion only if the qubit was in the down state at the 
time of measurement. Although this transition is forbid- 
den at the electric dipole level because of the spin fiip, 
spin-orbit coupling allows for a nonzero transition rate, 
as described below. 

Crucially, this architecture also allows for rapid ini- 
tialization of the qubit. Consider exposing a random 
ensemble of qubits to radiation of frequency 1^12. An 
electron in the (l,i) state will be excited to the level 
(2, t) and experience a relatively fast relaxation to the 
ground state (1, t), compared to a spin-fiip relaxation to 
the level (1, i). The net result is to rapidly polarize and 
thereby initialize the qubits. By varying the gate voltages 
(and thus 1^12) on individual dots we ensure that only de- 
sired qubits are brought into resonance. Clearly, under- 
standing the competition between the three time scales, 
l/jz/j, l/F/, and Ti, as a function of material parameters 
and gate potentials is the key to utilizing this device for 
readout and initialization. Note that Ti represents a 
thermal initialization time from (1,1) to (1,T)- Since Ti 
is a decoherence time, we require <C Ti. The Rabi 
oscillation frequency vji depends on the incident inten- 
sity and is therefore controllable, within limits. Robust 
measurement requires that many Rabi oscillations occur 
before orbital decay: l/T/. 

In order to understand the competition of time scales, 
we introduce an rms interaction energy, time averaged 
(as indicated by {}) over an optical cycle, which ex- 
presses the strength of the interaction in the electric 
dipole approximation. This defines the Rabi oscilla- 
tion frequency \hvR\^ = {\V^^\'^}. Here, V^^ ^ 
{—ehEo/m* Ei2)£ ■ P is the dipole term in the Hamilto- 
nian, |-BoP is twice the mean value of |E(t)p averaged in 
time, and e is the polarization unit vector. The electric 



field Eq inside the semiconductor with dielectric constant 
e = SrSo is related to the intensity of the external radia- 
tion / by i?o = Y^2//ceov^. 

The dipole Hamiltonian does not fiip the spin directly, 
but spin-orbit coupling causes each qubit state to be 
an additive mixture of up and down spin. In the 2D 
limit, the spin-orbit (so) Hamiltonian is dominated by 
the bulk [Dresselhaus (D)] and structural [Rashba (R)] 
inversion asymmetry terms, iJgo Hd + Hr, where 
Hd = P{Py(^y - PxCTx) and Hr = a{pxO-y - PyCTx)- In- 
cluding Hso perturbatively gives a nonzero dipole matrix 
element, and for light polarized in the y direction, the 
readout frequency is given by 

/-^-a^|(2|ya,|i)| 
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Note that the applied radiation need not be circularly 
polarized for this readout scheme. The Dresselhaus and 
Rashba parameters a and /? depend on intrinsic mate- 
rial properties, device design, and external electric field. 
For GaAs, both theoretical and experimental values vary 
widely O: a = 1 - 1000 m/s and (i = 1000 - 3000 
m/s. In a centrosymmetric crystal like silicon which has 
no bulk inversion asymmetry, /3 = 0. The one known 
data point for a SiGe two-dimensional electron gas gives 
a ~ 8 m/s 17], which we use in our estimates below. 

The relaxation of our quantum dot to its ground state 
enables spin polarization, but limits or even inhibits read- 
out if it occurs too quickly. This problem has been ad- 
dressed by Khaetskii and Nazarov in GaAs quantum dots 
[1^. In silicon, where there is no piezoelectric interaction, 
we calculate the relaxation rate via the Golden Rule with 
the usual deformation potential electron-phonon Hamil- 
tonian At sufficiently low temperatures (T < \K\ 
optical polar phonons and multiphonon processes do not 
contribute. By considering only longitudinal phonons, 
with dispersion lo = wjq, and using the long-wavelength 
approximation e*'' '" ~ 1 -f ik • r, we obtain the orbital 
decay rate due to electron-lattice coupling: 
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where p is the mass density, and "E.^ and S„ are the 
deformation constants. In strained systems, transverse 
phonons can also be important, as we shall discuss in a 
future publication [2o| . 

We perform a numerical analysis to obtain perfor- 
mance characteristics for the proposed measurement 
scheme. The numerical techni ques used are an exten- 
sion of those used in Refs. [sl lllj. The gate poten- 
tials, the electronic orbitals, and their corresponding im- 
age potentials (arising predominantly from the metallic 
gates) are computed self-consistently by a combination 
of three-dimensional finite element and diagonalization 
techniques. Specifically, we determine the readout oscil- 
lation frequency Eq. (Q, the orbital decay rate Eq. ||2Jl, 
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FIG. 2; The readout device studied in Figs. 3 and 4. Many 
similar devices were modeled in this work. 

and the coupling sensitivity of the qubit electron to an 
integrated SET. We have modeled numerous devices and 
a full account will be published elsewhere |2^. One of 
the most promising is shown in Fig. |3 This is a 6 nm 
strained silicon quantum well sandwiched between layers 
of strain-relaxed silicon- germanium (SigsGcis). The bot- 
tom barrier (30 nm) separates the quantum well from a 
grounded metallic back-gate. The top barrier (30 nm) 
separates the quantum well from lithographically pat- 
terned Schottky top-gates, whose voltages can be con- 
trolled independently. We consider negative gate poten- 
tials, which provide lateral confinement of the quantum 
dot through electrostatic repulsion. The main features of 
the quantum dot are that it is narrow, long, and slightly 
asymmetric. The narrow feature ensures that the ex- 
cited orbitals are non-degenerate, so a microwave field 
with narrow line width will not induce unwanted transi- 
tions. The long feature ensures that the energy splitting 
Ei2 (and thus F/) will be small. The asymmetry ensures 
that excited orbitals will provide charge motion. Figure 3 
shows the confinement potential for the device of Fig. 2; 
the asymmetry of the dot potential is apparent. 

Our most interesting results are obtained in the regime 
where the gate and image potentials are comparable in 
size. The effective confinement potential, including im- 
ages, is rather complicated due to the inhomogeneous 
gate arrangement. Indeed, it is a very poor approxima- 
tion to neglect screening in this system. The computed 
wavefunctions for our dot are shown in Fig. 3. 

We have also incorporated a SET into our proposed 
readout device, as shown in Fig.|21 The island of the SET 
is sheathed by a thin, 2 nm layer of silicon dioxide, and 
tunnel-coupled to adjoining source and drain gates. A 
third, capacitively-coupled gate is placed nearby, to pro- 
vide full control over both the charge and the potential of 
the island. When the device is operated in the Coulomb 
blockade regime, it becomes a sensitive electrometer [2l| . 
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FIG. 3: Electrostatic confinement potential and qubit elec- 
tron wavefunctions. The bare potential shown here is ob- 
tained in the quantum well along the device symmetry line, 
a; = 0. The contour plots show the electron probability den- 
sities in the x-y plane. 



In the gate arrangement of Fig. [3 the SET plays two 
roles. First, voltage control of the island allows us to 
vary the size of the dot and thus the energy splitting 
Ei2. Second, capacitive coupling to the quantum dot en- 
ables detection of its orbital state. The scheme works as 
as follows. Since the dot orbitals are spatially distinct, 
they induce different amounts of charge on the SET is- 
land. Consequently, transport currents through the SET 
will reflect the orbital states of the dot. The device ex- 
hibits optimal sensitivity if biased at the half-maximum 
of the conductance peak. The third SET gate has been 
introduced specifically to adjust this working point. As 
expected for our geometry, we find that the SET couples 
most strongly to the excited electronic orbital. 

For the device shown in Fig. [5] we obtain an energy 
splitting of £^12 = 0.129 meV = 31.2 GHz between the 
two lowest orbital states, and the dominant matrix el- 
ements \{l\y\2)\ = 48 nm and \{2\ydy\l)\ = 3.6. From 
these results we obtain the readout oscillation frequency 
lyn = 5.5 X 10^\/l Hz (for microwave intensity / in 
units of W/m^), and the orbital decay rate for spon- 
taneous phonon emission F/ = 12.7 MHz. The orbital 
decay rate for emission of a photon is very much less: 
F^ = e^Ef2\{l\y\2)\^V^/STrh'^c^el^^ = 6.5 Hz, and there- 
fore is not a limiting factor in this scheme. All results 
are obtained at an ambient temperature of T = 100 mK 
and a magnetic field of H = 0.05 T = 1.40 GHz. The 
dielectric constant of Si is £ = 11.9eo- Note that the 
natural widths of the states are small compared to their 
separation. The charge on the SET island was com- 
puted by integrating n • D over the surface. The excess 
or induced charge of the excited quantum dot orbital, rel- 
ative to its ground state, was found to be AQ = 0.052e. 
The corresponding center-of-charge motion in the dot is 
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FIG. 4; Summary of operation parameters for the device 
shown in Fig. 2. The shaded region shows the operating 
range, with constraints set by the sensitivity of the SET and 
the spontaneous decay rate of the upper readout state. The 
dotted Hne shows the actual operating range of the device 
when the intensity of the incoming radiation is left variable. 
A likely working point of pr, marked by x, corresponds to an 
intensity I = 7.2 nW//im^. 



Ay — 4.3 nm. However, the charge motion does not 
track closely with AQ; the latter is determined primarily 
by the capacitive coupling between the dot and the SET. 
Finally, we find that by changing the bias voltage on the 
SET by 20% (thus reducing the dot size), the excitation 
resonance frequency vi2 changes by 8 GHz. 

The prospects for initialization and readout in this 
scheme are partially summarized in Fig. 4. Theoreti- 
cal considerations of shot noise [21I place an upper 
bound of about 4 x 10^^e/\/Hz on the detection sensi- 
tivity for charge induced on the island of an optimized 
rf-SET, as a function of the measurement bandwidth. 
Similarly, the decay rate of the excited electronic orbital 
places a lower bound on the readout oscillation frequency 
Vfi. The latter is a function of microwave intensity. By 
directing 70 pW of microwave power onto a dot of size 
~ 0.1 X 0.1 ^iTc? , as consistent with low-temperature 
transport experiments, we would obtain the working 
point marked in Fig. 4, with l/vji ^ l/F/. Sample heat- 
ing outside the dot can be minimized by using integrated 
on-chip antennas to focus the power 23]. We must also 
satisfy the requirement Tg ^ l/E/ = 78.7 ns where Tg 
is the spin decoherence time. Theoretical estimates for 
Ti in Si quantum dots exceed 1 ms Experiments 
on dots have yet to be performed, but spins of electrons 
in impurity donor states should have similar lifetimes. 
Recent work on electron spins in ^*Si:P [loj would im- 
ply that Ts > 60 ms. Hence the crucial hierarchy of 
timescales I/i'r ^ 1/^/ <C Ts should be achievable. 



In summary, we have proposed, simulated, and an- 
alyzed a scheme for measuring single electron spins in 
a quantum dot via spin-charge transduction. Adequate 
charge detection sensitivity is provided by an integrated 
rf-SET. The scheme also enables fast initialization. 
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